Bose-Hubbard model on a star lattice 
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We analyze the Bose-Hubbard model of hardcore bosons with nearest neighbor hopping and 
repulsive interactions on a star lattice using both quantum Monte Carlo simulation and dual vortex 
theory. We obtain the phase diagram of this model as a function of the chemical potential and the 
relative strength of hopping and interaction. In the strong interaction regime, we find that the Mott 
phases of the model at 1/2 and 1/3 fillings, in contrast to their counterparts on square, triangular, 
and Kagome lattices, are either translationally invariant resonant valence bond (RVB) phases with 
no density-wave order or have coexisting density-wave and RVB orders. We also find that upon 
increasing the relative strength of hopping and interaction, the translationally invariant Mott states 
undergo direct second order superfluid-insulator quantum phase transitions. We compute the critical 
exponents for these transitions and argue using the dual vortex picture that the transitions, when 
approached through the tip of the Mott lobe, belong to the inverted XY universality class. 

PACS numbers: 75.10.Jm, 05.30.Jp, 71.27.+a, 75.40.Mg 



I. INTRODUCTION 



The study of microscopic models which may lead to 
exotic quantum phases has been carried on for a long 
time in condensed matter physics. Recently, one such 
model system, namely, the two-dimensional (2D) Bose- 
Hubbard model on a lattice, has received a great deal 
of attention. One of the reasons for this renewed atten- 
tion is the possibility of experimental realization of such 
a model using cold atoms trapped in optical lattices. ^'^ 
However, the Bose-Hubbard model is theoretically inter- 
esting in its own right. In particular, it has recently been 
pointed out that the superfluid-insulator transitions for 
fractional Boson flUing factors in this model may be of 
non Landau-Ginzburg type in the sense that the low en- 
ergy theory for these transitions cannot be described in 
terms of order parameter fields of the phases in either 
side of the transitions.'^ Instead, as pointed out in several 
studies of the modelj^i^i^ the transition is aptly described 
in terms of vortices which are topological excitations of 
the superfluid and whose condensation ultimately leads 
to destabilization of the superfluid phase in favor of insu- 
lating Mott phasesi^ Such a dual vortex theory provides a 
list of possible competing Mott phases via general sym- 
metry requirements of the underlying lattice.'^ In par- 
ticular the geometric frustration induced by the lattice 
structure, which plays a key role in determining the na- 
ture of these competing Mott phases, is quite naturally 
described by the dual vortex theory. 

Another, more direct and quantitative, approach 
to studying these Bose-Hubbard models on a lat- 
tice has been numerical quantum Monte Carlo (QMC) 
simulationsi^iii^ These QMC studies provide us with 
quantitatively accurate phase diagrams of the model. 
Further, they also furnish direct information about both 
the nature of the phase transition and the correlation 



functions in the Mott phasCiSjl'^ These features, along 
with the possibility of accessing much larger system sizes 
than possible in exact diagonalization studies, make them 
the numerical method of choice for the lattice Bose- 
Hubbard models. A combination of the dual vortex the- 
ory and QMC simulation has been recently used to study 
the phases of Bose Hubbard model on triangular'''* and 
Kagome latticesj^i^ As noted in these works, the Bose- 
Hubbard model can also be mapped onto a spin-1/2 XXZ 
modcl^ leading to interpretation of the obtained results 
in terms of both bosons and quantum spins. 

More recently, there have been several studies of 
the antiferromagnetic Heisenberg model on a star 
latticeJ^iiiii^ This lattice, shown in Fig. [U can be con- 
structed by expanding the sites of a hexagonal lattice 
into triangles. It consists of two topologically inequiva- 
lent bonds, triangular and expanded, as shown in Fig. [TJ 
It has been argued in Ref. [Ill, using exact diagonalization 
studies, that the ground state of the Heisenberg model 
on this lattice is a paramagnetic valence bond crystal 
(VBC). Such a VBC state is characterized by enhanced 
antiferromagnetic correlations along the expanded bonds 
indicating singlet formation along those bonds. Such 
studies have also been extended for the anisotropic anti- 
ferromagnetic Heisenberg model where the exchange cou- 
pling j"^ on triangles is different from the coupling on 
expanded links.'*' It has been argued in Ref. |lQ that there 
may be another VBC state for j"^ > 1.3J^ which consists 
of 18-site star pattern. Possible spin liquid and valence 
bond crystal phases are also studied in an Sp{N) gen- 
eralized model. '^ These studies, till date, have not been 
extended to other spin models such as the XXZ model. 

In this work, we study the Bose-Hubbard model, or 
equivalently, the XXZ model on a star lattice. The 
Hamiltonian of the model can be represented in terms 
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of hardcore bosons as 
Hb = - [t''h\b, + h.c) + y"'n,n, - M 

(1) 

where t'^^{V'^^) is the hopping amphtude (interaction 
strength) of the Bosons between sites i and j, bi denotes 
annihilation operator of the boson at site i, rii = blbi is 
the number density operator for the Bosons at site i, (ij) 
indicates that sites i and j are nearest neighbors of each 
other, and /i is the chemical potential. This model, in 
contrast to the Heisenberg model on frustrated lattices, 
is amenable to QMC studies. In what follows, we shall 
allow for different amplitudes of hopping amplitude and 
interaction strengths: f^^ — and V'^^ — for ex- 
panded bonds and f^^ = and V^^ — for triangular 
bonds. In this work, we shaU set /t^ = /V'^ . We 
note at the outset that this model can be mapped onto an 
XXZ model in a magnetic field via a Holstein-Primakoff 
transformation^ and yields 

i/xxz = -\Y.jr\sts-+s-s^] 

iv) 

+ YJl'StS^^+KY.^l (2) 

where jf = V-^ , Jl^ = V^^ , h^^ {^i- 1/2), and the spin 
operators are expressed in terms of the boson operators 
as 5'+ = bl,S- = bi,Sf = b\bi - 1/2 in the leading 
order in 1/S. We shall use these spin and the boson 
representations of the model interchangeably throughout 
the paper. Note that half filling in boson language means 
zero magnetic field in spin language. 

The central results reported in this work are the fol- 
lowing. First, using QMC simulation for sufficiently 
large systems {L < 60) and low temperatures (/3~^ < 
O.OOIJ^), we obtain a phase diagram of the Bose- 
Hubbard (XXZ) model as a function of the fi/V^ {h^/ J^) 
and t^/y^(jf /Jf ) for fixed ratios /V^{t^ /t^). We 
find that there are three distinct Mott phases at bo- 
son fillings 1/2 and 1/3. Second, using QMC simula- 
tions, we compute the equal-time spin-spin correlation 
functions and the real space bond-bond correlation func- 
tions for the bosons (spins). From these studies, we 
demonstrate that two of these Mott phases, which oc- 
cur for = at 1/2 and 1/3 fillings, are transla- 
tionally invariant and do not exhibit density-wave order 
(magnetization). These Mott states, in stark contrast 
to their counterparts in square, triangular or Kagome 
lattices^iSii^ are found to have resonating valence bonds 
(RVB) along either the triangle or the expanded links. 
The third Mott phase has coexisting density-wave (Neel) 
and RVB orders and occur for ^ and at 1/2 fill- 
ing. These QMC results regarding the nature of the Mott 
phases are also supported by qualitative symmetry-based 
analysis using dual vortex theory. Third, the transla- 
tionally invariant Mott states occuring for = are 



found to undergo a second-order superfluid-Mott insula- 
tor quantum phase transition with increasing m 
the spin language, this transition corresponds to a shift 
from Sx to Sz ordering if Sz ordering is present or from Sx 
ordering to a paramagnet if such an ordering is absent). 
This is also in contrast to analogous studies on square, 
triangular and Kagome lattices^i^i^ where QMC simula- 
tions found evidence of either a direct first-order transi- 
tion between the superfluid and the Mott phases or an 
intermediate supersolid phase. A qualitative symmetry- 
based analysis using the dual vortex theory finds this 
transition to be in the inverted XY universality class, 
when the transition point is approached through the tip 
of the Mott lobes. This observation is also supported by 
a finite-size scaling studies using QMC simulations which 
yields the dynamical critical exponent z and the correla- 
tion length exponent ly for the transition. 




FIG. 1: The star lattice. The 'triangle' bonds are denoted 
by thin lines and the 'expanded' bonds are denoted by thick 
lines. 

The organization for the rest of this work is as follows. 
In Sec. nil we elucidate our QMC results. We deal with 
the case Jj = Jf in Sec. Ill Al and follow it up with 
the study of the case Jj 7^ jf for 1/2 filling in Sec. 
IIIBI These numerical results are then compared with 
the analytical predictions of the dual vortex theory in 
Sec. mil This is followed by conclusion in Sec. |IVl 



II. QUANTUM MONTE CARLO 

In this section, we shall analyze the model using quan- 
tum Monte Carlo. We use a multi-site generalizationi^ 
of the stochastic series expansion (SSE) methodJ^ Here 
the basic lattice unit is a site and all its neighbors. Sim- 
ulations are performed for systems of linear size L = 
12, 18, 24, 36, 48, 60 with N = QL"^ sites at different tem- 
peratures (the lowest temperature is f3~^ — Jj_/1200). In 
Sec. Ill Al we study the isotropic case with equal hopping 
and interaction strength on the triangle and expanded 
bonds while the anisotropic case is studied in Sec. Ill BI 
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A. J'^ = J== 

1. Phase diagram 

The phase diagram obtained for = is shown in 
Fig. [21 There are three phases: the XY ferromagnet and 
two valence bond crystals denoted by VBCl and VBC2. 
In the boson language, they correspond to the super- 
fluid and Mott phases respectively. The characteristics 
of these Mott phases and their transition to the super- 
fluid phase is discussed in Sec III A 2] and ULA 31 in details. 
We note here that the Monte Carlo scans are performed 
only along three lines shown in the phase diagram so that 
the phase boundaries are approximate. 
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FIG. 3: Scaling of the superfluid density ps ior z = 1, (i = 
Hal} filhng 16L/Jx and h^/J± = 0. Lines guide the eye. In this and all 

other figures, error bars are smaller than the symbol size if 
not visible. 





AKL 



FIG. 2: The schematic — phase diagram from Monte 
Carlo simulations. The phase boundaries are denoted by thick 
solid lines. Monte Carlo scans are denoted by dashed lines. 



FIG. 4: Data collapse of the superfluid density ps for z = 1, 
V = 0.67, (J^/Jx)c = 3.0783, /3 = 16L/Jx, and /i^/Jx = 0. 



In this subsection, we focus on the model in the ab- 
sence of an applied magnetic field, i.e. at 1/2 filling for 
bosons. There is a continuous transition from the super- 
fluid phase to an insulating phase at Jz/ Ji_ ~ 3.0783. In 
the vicinity of a continuous transition, the spin-stiffness 
of the XY ferromagnet (or superfluid density in the bo- 
son language), ps, which is measured through winding 
number fluctuations, scales as 

p,^L-^Fp^{L'/''iK,-K),(3/L^), (3) 

where Fp^ is the scaling function, L is the linear system 
size, z is the dynamical critical exponent, v is the correla- 
tion length exponent, SKc = Kc — K = {Jz/ J±)c — Jz/ J± 
is the distance to the critical point, f] — (ksT)^^ is the 
inverse temperature, and fcs is the Boltzmann constant. 
It follows from the above finite size scaling relation (Eq. 
[3]) that the curves for different systems sizes should cross 



at the transition point when psL^ is plotted as a func- 
tion of the coupling constant for fixed (or for large 
enough /? to ensure the ground state convergence) . It also 
follows from Eq. [3] that the curves for different system 
sizes should collapse onto a universal curve for appropri- 
ate values of v and {Jz/ J±)c when pgL^ is plotted as a 
function of SKcL^^'^ ■ The data scale well with the dy- 
namical critical exponent z — \. In Fig. [31 psL is shown 
as a function of the coupling constant. The curves for 
different system sizes cross at a distinct point. The data 
collapse is shown in Fig. [Hand leads to a critical exponent 
V = 0.67. 

To address the nature of the insulating phase for Jz S> 
Jj^, we study the equal time and static spin structure 
factors that are given by 

5(q) = L^Sl^S^r), xicD^L'Q dr5t,5qo^(4) 
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Thus the real-space bond-bond correlation function also 
confirms the formation of singlets along the expanded 
bonds. This lead us to conclude that the most proba- 
ble candidate for this insulating phase is VBC with no 
symmetry breaking and is analogous to the VBC state 
described in Ref . VLVi We denote this phase as VBCl. 



3. 1/3 fillmg 



FIG. 5: The equal time spin-spin structure factor (left panel) 




FIG. 7: Scaling of the superfluid density ps for z — 2, f3 = 
FIG. 6: The correlation function Cb(ri - r^) between the L^/3J±, and J^/J± = 4. Lines guide the eye. 
bond indicated by a large text and the other bonds for L = 24, 
J^/Jj^ = 4, h^/J± = 0, and ksT = O.Of Jx. 



where S^^t = ( 1 /-^^ ) X]j '^'fr 6xp(«q ■ r^). The structure 
factors in the insulating phase are shown in Fig.O There 
are no sharp peaks which is a clear indication of ab- 
sence of magnetic order. The short-range bow-tie fea- 
tures, found to be present in the structure factor, are 
likely to be a remnant of classical dipolar correlations at 
finite temperatures, that are known to arise in classical 
Heisenberg models on various frustrated lattice. 

The nearest neighbor antiferromagnetic spin-spin cor- 
relations are strongly enhanced along the expanded 
bonds (bonds that connect triangles) signaling singlet 
formation along those bonds: {[S^S%) w 2.92([S'^S'^]t) 
for J,/J± = 4, hJJ± = 0, and T = O.OIJ^. To ver- 
ify this, we further compute the real space bond-bond 
correlation function that is given by 



Cb(r^ - r^-) = J B^rdr J Bsrdr 



(5) 



where BQ-(ij),i- = J±{SfSj + SfS^) is the off-diagonal 
bond operator (at imaginary time r) of the bond a con- 
necting spins i and j. A plot of Cb{r^ — rs) in Fig.[B]shows 
that the off-diagonal bond operators are distributed uni- 
formly on the expanded and triangle bonds with the ma- 
jority of operators on the expanded bonds with a ratio 
10.1/2.95 « 3.42. There is no any other bond order. 
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FIG. 8: Data collapse of the superfluid density ps for z = 2, 
V = 0.5, (/!.^/Jx)c = 2.37568, P = /iJ^. AK = 
{hz/J±)c ~ hz/ Ji_ at fixed Jz/Ji_ = 4. 

In this subsection, we focus on the phase diagram away 
from the 1/2 filling, or, in the spin language, in the pres- 
ence of an applied longitudinal magnetic field. As shown 
in Fig. [2l there is an additional VBC phase (denoted as 
VBC2) at the filling of 1/3 (spin magnetization equals 
— 1/6). The transition from the XY ferromagnet phase 
to the VBC2 phase is continuous. Within the system 
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FIG. 9; The equal time spin-spin structure factor (left panel) 
and the static structure factor (right panel) for L = 24, 
Jz/J± = 4, hz/J± — 3 and ksT = O.OOSJx. Ferromagnetic 
peaks due to the uniform background magnetization are sub- 
tracted. The axes range from —Att to Att. 




the triangle bonds in the VBC2 phase forming trimers. 
There is no other bond order. We can conclude that the 
VBC2 phase is a quantum paramagnetic phase without 
any symmetry breaking. 



B. 



J'^ / at /i^ = 




FIG. 11: The schematic ^ phase diagram at 1/2 filling 
from Monte Carlo simulations. The phase boundaries are 
denoted by thick solid lines. 



FIG. 10: The correlation function Cb(ri — rs) between the 
bond indicated by a large text and the other bonds for L = 24, 
Jz/Jj^ = 4, hz/J± = 3, and ksT = 0.01 Jx. 



size and temperatures that we have studied, we have not 
found any signatures of a first order transition such as 
double peaked histograms. As shown in Fig.[7]and Fig. [51 
the data scale reasonably well with z = 2 and ly = 0.5. 
It is worth mentioning that the scaling is not as good 
as that for the transition to the VBCl phase described 
in the previous section. The quality of the scaling is 
probably limited by the fact one needs data of very high 
accuracy at extremely low temperatures and very large 
system sizes (in order to reach the scaling regime given 
by Eq. [3]) and that is beyond our computational facilities 
at the present time. 

As shown in Fig. [9l the VBC2 phase does not exhibit 
magnetic order. The connected spin correlation func- 
tion is shorter-ranged than the spin correlation func- 
tion in the VBCl phase. This can be deduced from 
the larger bow tie width in Fig. [9] compared to Fig. [5l 
The nearest neighbor antiferromagnetic spin-spin corre- 
lations are slightly weaker along the expanded bonds: 
{[S'S%) « 0.78{[S'S']t) for J,/J^ = 4, h,/J^ = 3, and 
ksT = 0.01J_L. As shown in Fig.fTOl the ratio of the num- 
ber of the off-diagonal operators on the expanded bonds 
to that on the triangle bonds is 0.45/1.71 « 0.26. Thus, 
in contrast to the VBCl phase, the spins resonate along 



The schematic J'^ ^ phase diagram at 1/2 filling 
is shown in Fig. [TT] There are again three phases: the 
XY ferromagnet (superfluid in the boson language), the 
valence bond crystal phase from Sec. Ill Al (VBCl). and a 
phase that is characterized by both VBC order and Neel- 
like magnetic order and is denoted as VBC-Neel. The 
latter two phases are Mott states in the boson language. 
Monte Carlo scans are performed only along a few lines 
so that the phase boundaries are approximate. We have 
not attempted to determine the nature of phase transi- 
tions. However, as shown in Fig. 1121 we find a narrow 
region with finite superfluid density between the VBCl 
and VBC-Neel phases even for large Jj/Jj. 

In Fig. [131 we show the real space bond-bond corre- 
lation function for different values of Jl'/jJ. There are 
singlets on the expanded bonds for large Jf/Jj (VBCl 
phase) and there are resonating triangles (trimers) for 
smah /Jj (VBC-Neel phase). The VBC-Neel phase 
also exhibits long range magnetic order. In Fig. 1141 
we show the finite size scaling of the equal time struc- 
ture factor given by Eq. [31 at the ordering wave vector 
Q = (27r, 27r). Note that the structure factor vanishes at 
Q = (0, 0) due to geometrical factors. The structure fac- 
tor divided by the number of lattice sites clearly scales to 
a finite value in the thermodynamic limit indicating long- 
ranged magnetic order. The triangles of the star lattice 
form the bipartite hexagonal lattice. The structure of 
the real space correlations is such that the spins on two 
different sublattices of the hexagonal lattice (belonging 
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FIG. 12: The superfluid density ps as a function of Jf/Jj 
for Jj/Jj = 10, L = 12, and ksT = O.OUl. 
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FIG. 13: The correlation function Cb(ri — r^) between the 
bond indicated by a large text and the other bonds in the 
VBCl phase at Jf /Jj = 0.6 (upper panel) and in the VBC- 
Neel phase at Jf/Jj = 0.1 (lower panel) for Jj/Jj ~ 10, 
L = 12, and ksT = O.OljJ. 



to up and down triangles in Fig. [T]) have antiferromag- 
netic correlations as in a Neel-like ordered state. This 
long-ranged magnetic order corresponds to density-wave 
order in the boson language. 

To understand the VBC-Neel state in more detail, con- 
sider isolated triangles (J^ = 0). There are two degen- 
erate ground states per triangle. These are resonating 
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FIG. 14: The equal time structure factor at the ordering wave 
vector Q as a function of the inverse system size for jj/ jj = 
10, L = 12, and ksT = O.OUl. 



trimers with the total spin 5^ — —1/2 
1 



IV'i 



V3 



(I- 



and the total spin = 1/2 



1 



where -|-(— ) indicate sites with — l/2(— 1/2). In the 
boson language, and IV'2) correspond to one and two 
bosons per triangle respectively. At 1 /2 filling, the num- 
ber of spin up triangles is equal to the number of spin 
down triangles. For = 0, those triangles can be ar- 
ranged arbitrary on the star lattice and the ground state 
is extensively degenerate. However, finite selects a 
Neel state with respect to the total spins of the triangles 
because the triangles form a bipartite lattice and the ef- 
fective interaction between them is antifcrromagnetic. A 
similar state is also found in the dual vortex theory anal- 
ysis, see the next section. 



III. DUAL VORTEX THEORY 

In this section, we shall obtain an analytical under- 
standing of the nature of the Mott phase and the quan- 
tum phase transitions from them to the superfluid phases. 
Throughout this section, we shall restrict ourselves to the 
isotropic case = . 

The derivation of a dual vortex action starting from 
the Bose- Hubbard model (Eq. [T|) has been elaborated in 
Refs. dQH. The vortices are described in terms bosonic 
field ipb and a dual gauge field Ab^ which lives on the 
sites b and links fj, of the dual lattice respectively. A 
duality analysis of the Bose- Hubbard model then leads to 
an effective dual action which can be expressed in terms 
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FIG. 15: The star lattice and its dual. Notice that the A sites 
are in the center of the hexagon of the dual lattice while the 
B and the C sites are in the center of the triangle. 



of the vortices and the gauge fields aa^ 



Z 

Sd 



J 'DA J Vecxp{-Sd) 

b 
b 



(6) 



where V'h are the vortex field living on the site b of the 
dual lattice, Ai,^ is the U(l) dual gauge field such that 
eri^xAijAbx = rii where nt is the physical boson den- 
sity at site i, J2p denotes sum over elementary plaquette 
of the dual lattice, denotes lattice derivative along 
H = x,y, T, and / is the average boson density. Here Uv 
is the vortex fugacity and r, u, and e denotes parame- 
ters of the dual action which can not be directly mapped 
onto those of Hi, since Sd is not self-dual to the boson 
action obtained from Hh. Therefore we cannot, in gen- 
eral, obtain a direct mapping between the parameters of 
the two actions, except for identifying the magnetic field 
seen by the vortices eruX^uAbx as the physical boson 
densityi^'^'^ In the remainder of the paper, we shall clas- 
sify the phases of this action based on symmetry consid- 
eration and within the saddle point approximation where 
the gauge fields are pinned to their saddle point values. 

The star-lattice and its dual is shown in Fig. [121 The 
first step towards understanding the phases of Sd within a 



saddle point approximation amounts to solving the Hopf- 
stadter problem for the vortices on this dual lattice shown 
in Fig. [m We note from Fig. [15] that the lattice dual to 
the star represents a dice lattice with all diagonals con- 
necting the A sites joined to each other. The Hamiltonian 
of the vortices on such a lattice is given by 



(y) a,l3=A,B,C 



h.c 



(7) 



where ipia = '<Pa{o-x, o-y) denotes the annihilation opera- 
tor for vortex fields at i = {ax, ay) and 7^- is the dual 
magnetic flux, which in the gauge A = H(0, x), is given 

by 



lij = 27r/(2a^/a + Xa){yj ~ y'i) 



(8) 



where we have used Xi — a^, Xj = a^ + AqO, y' = 
2j//(V3a), / — -ffV3a^/40o is the flux passing through 
an elementary triangle in units of basic flux quanta (/>o, 
and Aq is a number that depends on the sublattice index 
a — A,B, C. Here the dual flux passing through an el- 
ementary plaquette is double of that of the dice lattice. 
This can be inferred from the fact the present lattice has 
two sites of the star lattice in each rhombus (as opposed 
to one site of the Kagome lattice in each rhombus) of the 
dice latticei^ii^ 

The Schrodinger equation for the vortex flelds can 
be constructed from the Hamiltonian. With our choice 
of the gauge, we can write ipa{cix,o,y) — e"*"^ i^dax), 
where Ky — 2aky/\/3 S (0, tt) since wc are restricted 
within the first Brillouin Zone, and a — A, B,C' rep- 
resents the inequivalent sites of the dual lattice. We 
thus obtain, writing energy e in units of t and defining 
0± = 2'Kf{2ax/a ± 1/2) + Ky, 

- ei^siax) = ipA{ax + a) + 2ipA{ax - a /2) cos (j)+{9) 
-e^cia-x) = -tpAiax - a) + 2-11; A{ax + a/2) cos 4>i-10) 
-etpA^ax) = ^pB{ax + a)+ipc{ax-a) 
+2tl;B{ax - a/2) cos((/)_) 
+2ipciax + a/2) cos((?!)+) 
+2tpA{ax + 3a/2) cos(27r/ + (j)+) 
+2ipAiax - 3a/2) cos(-27r/ + (/)_) 
+2il;A{ax) cos{4)+ + (j>^) (11) 

Note that if we ignore the last three terms in the RHS 
of Eq. [TT] which involves ipA, we get back the dual Hop- 
stadter equation for the dice latticei^ii^ To solve for e, 
we substitute Eqs. M and [TU] in Eq. [TT] and get, using 
Ux = 3rna/2 for all A sites. 



6)V'A77 



(4cos(2^/)-2e)C™ 

IpArn+l COs(67r/(m + 1/2) -f Ky) 

+ijjAm.-i cos(67r/(m - 1/2) + Ky) 



+tpAm COs(127r/TO + 2Ky 



(12) 



We first consider the case / = 2/3, which is identical to 
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FIG. 16: Mott state at / = 2/3. The circles indicate a trian- 
gle with either two spin up and one spin down on the three 
vertices or a trimerized triangles indicted by circles when su- 
perposition of the spins/bosons due to quantum fluctuations 
are allowed. The state is translationally invariant and agrees 
with the VBC2 state predicted by QMC when quantum fluc- 
tuations are included. The corresponding ground state for 
f — 1/3 can be obtained by simply flipping the spins. 



the / = 1/3 filling considered in Sec. Ill A 31 Here Eq.[T2] 
reduces to 

(e2-6)V'Ak = -2(1 + e)Aki^Ak 

v4k = 2 cos(fc^) cos(K.y) -f cos(2kj,) (13) 

where we have taken Fourier transform with respect to 
m and k'^ = 3kxa/2 e (0,7r). This has the solution 



(^k - 1) 



1/2 



(14) 



so that the minima of the vortex spectrum occurs at 
{k'^,Ky) = (0, 0), (tt, tt) and corresponds to e = —6. Also, 
substituting the values of k'^ and Ky in Egs.l^and llOl we 
find that ^^(fc^) = = tpcik^) so that the eigenfunc- 
tions corresponding to {k'^,Ky) — (0, 0), (tt, tt) are given 

by 

^i(m, n) = (1, 0, 0) n) = (1, 0, 0)e'^(™+") (15) 

where we have used ay — \/3na/2 for the A sites. Note 
that since m -I- n is always even for A sites "01 = "4^2 and 
hence the theory has a single vortex field which should 
mediate the transition. Also note that the Mott state is 
expected to be uniform since is basically a constant. 
Thus the simplest state compatible with these require- 
ment at 2/3 filling is shown in Fig. [121 Here each trian- 
gle denoted by a circle, within mean-field, has two occu- 
pied and one empty sites leading to a net occupancy of 2 



bosons per triangle. Equivalently, in the spin language, 
this corresponds to two spin-up and one spin-down sites 
at every triangle leading to a net magnetization of 1/2. 
However, it is indeed possible that inclusion of quantum 
fluctuations will make the bonds within the triangle to 
resonate leading to the trimerized VBC2 state obtained 
in QMC studies. 

The theory of transition pertains to a single vortex 
field in the presence of a fluctuating dual gauge fleld 
and is thus belongs to the inverted XY universality class 
which has z = 1 and f = 2/3^^ These exponents are 
the same as their counterparts for models in the 3D XY 
universality class. However, the key difference is that 
this exponent is obtained for a fully interacting model 
in = 2 -|- 1 in the strong interacting regime. This 
is in contrast to D = 3 -I- 1 dimensional systems where 
one expects the transition to be fluctuation-driven first 
order. We note that this expectation, which was ini- 
tially derived using an e = 4 — d expansion method, is 
not valid for D = 2 + 1 quantum systems where such 
transitions remain continuousJ^ Note that the quantum 
phase transition described by this dual vortex theory re- 
quires a fixed density across the transition and hence is 
valid when the transition is approached via the tip of the 
Mott lobe. The QMC study of Sec. Ill A 31 approaches 
the transition from the side of the lobe and hence gets a 
different z. 

Next we come to case of / — 1/2. Substituting / = 
1/2, in Eq.[Il one gets 



(e^ - 6)lpAm. = 
Cm — 



-4(1 + e/2)Cm 
-■ipAm+i sin(37rTO + Ky) 
+ipAm-i sin(37rm + Ky) 

+i)Am C0s(2Ky) 



(16) 



Thus here we need to distinguish between the sites which 
have m as even and odd integers. Denoting the corre- 
sponding fields as -0*^ and respectively, we find that 

[e^- 6 + 4(1 + e/2)cos(2K^)]i^^k + Ck0Ak = 

[^-(;> + A{l + e/2)cos{2Ky)]rA^. + C^rA^. = 

Ck = 8isin(fc;)sin(Ky)(l-t-e/2) (17) 

Note that here k'^,Ky G (0,7r/2) since the periodicity in 
real space has been doubled. From Eq. [T71 we find that 

2 where 



e = -Ak±±V(^k±-2 



Ak± = cos(2Ky)±2|sin(fc^)||sin(Kj^)| (18) 

Thus the minima of the spectrum occurs at (k'^,Ky) = 
(7r/2,7r/6) with e = —3. Substituting the value of e in 
Eq. [13 we find — ii/ja- Also, substituting the values 
of iPa/4'a — h k'^ — ''■/2 and Ky = 7r/6, we find ipB and 
ipc from Eq. 1161 Finally, this yields the wavefunction 



^1 



(o) 



[l(-i),0,c(-CT)]e'""ie*""^/^ 



(19) 



where c = 2/3 , the coordinates of odd and even sites are 
taken to be (01,02) = {2m + l,2n -I- 1) and (01,02) = 



FIG. 17: A possible Mott state with all the connecting bonds 
of the triangles forming dimers. The dimerized bonds are 
shown by thick lines. This state can occur for n = 0, 1, 2 or 3 
and is identical to the VBCl phase predicted by QMC when 
quantum fluctuations are included. See text for details. 

(2m, 2n) and we haven't renormalized the wavefunction. 
Thus, the theory of transition again pertains to a theory 
of single vortex field in the presence of fluctuating dual 
gauge field and belongs to the inverted XY universality 
class with z = 1 and v — 2/3. This is compatible with 
the exponents obtained by QMC study in Sec. Ill A"2l 

Finally, we consider the possible Mott states for / — 
1/2. These are obtained from qualitative argument out- 
lined below and are shown in Figs. [171 andfTSl From the 
vortex wave function (Eq. [T9|l , we find that all the B sites 
of the dual lattice are equivalent. Thus the triangles of 
the star lattice whose centers are occupied by B sites of 
the dual lattice must have the same filling and can be 
filled with n = 0, 1, 2 or 3 bosons (or n up and 3 — n 
down spins per triangle in the spin language). This leaves 
3 — n bosons (or 3 — n up and n down spins) to be dis- 
tributed over triangles which contain a C site of the dual 
lattice. Such a distribution must have the requisite peri- 
odicity of the wavefunction i.e. a 4 x 12 unit cell which 
involves 4 A sites (4 hexagons of the real lattice) in the x 
direction and 12 A sites (12 hexagons of the real lattice) 
in the y direction. Thus we find a multitude of energy- 
equivalent mean-field states with 4 x 12 units cells which 
corresponds to different ways of filling these sites keeping 
the total boson density (or magnetization) fixed to 3 — n 
on the C sites. Similar to the case of XXZ model on a 
Kagome lattice analyzed in Ref. these states can su- 
perpose in the presence of quantum fluctuations leading 
to a pattern having n Bosons in every triangle which has 
a B site of the dual lattice at its center and 3 — n bosons 



FIG. 18: Another possible state at / = 1/2 for n = 2. The 
large (red) circles indicate B sites and the corresponding tri- 
angles have one dimerized bond and an up spin while the small 
(blue) circles mean a triangle with one dimerized bond and a 
down spin. A similar state for n = 1 can occur and can be 
obtained by simply interchanging the sizes (colors) of the B 
and the C sites. This state is identical to the VBC-Neel phase 
predicted by QMC when quantum fluctuations are included. 
See text for details. 

(on the average) in every triangles which has a C site of 
the dual lattice at its center. Now if we allow the spins 
on the joining links of these triangles to hybridize, we get 
a translationally invariant dimerized state, shown in Fig. 
[T71 where each of the connecting bonds of the triangle 
can form a singlet dimer. This state is analogous to the 
VBCl state obtained in QMC study. Another possible 
state, shown in Fig. [15] which correspond to n = 2 or 
n = 1, where the one bond in each of these triangles can 
hybridize (or form a valence bond), leads to the transla- 
tion symmetry broken state shown in Fig. [181 This state 
is analogous to the VBC-Neel state found in QMC stud- 
ies for J-^ 7^ . There may be other possible states and 
a full classification of all of them seems to be difficult. 
We point out that the classification of these Mott states 
necessarily requires incorporation of quantum-fluctuation 
induced superposition between possible mean-field states. 



IV. CONCLUSION 

In conclusion, we have presented a study of Bose Hub- 
bard or equivalently spin 1/2 XXZ model on a star lat- 
tice using both QMC and dual vortex theory. We have 
shown that for = , the model supports translation- 
ally invariant RVB Mott phases at / = 1/2 and / = 1/3 
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and have pointed out that these phases are different from 
their counterparts with broken translational symmetry 
in square, triangular and Kagome lattices. We have 
also shown that these phases, upon increasing the ra- 
tio of nearest neighbor hopping amplitude to interaction 
strength, undergo a direct second order quantum phase 
transition to a superfluid phase. We have identified the 
exponents of this transitions and shown that they belong 
to the (2-|-l)D inverted XY universality class with z = 1 
and V = IjZ when approached through the tip of the 
Mott lobe. When the transition is approached from the 
side of the Mott lobe for / — 1/3, QMC finds a second or- 
der transition with z — 2. Such clear signatures of second 
order quantum phase transitions is in contrast with the 
behavior of the model on square, triangular or Kagome 
lattice, where these transitions are either first-order or 



are accompanied by intermediate supersolid phases. We 
have also provided a phase diagram for the system at 1/2 
filling for 7^ and have demonstrated the existence 
of a Mott phase with coexisting density-wave (Neel) and 
RVB orders. 
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